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Inelastic (Raman) light scattering intensities for a 42-electron quantum dot under off-resonance 
conditions and in different spin and angular momentum channels are computed in order to test 
whether final collective states become the dominant peaks in the process. The results of the calcu- 
lations set a limit to the spread use of the off-resonant approximation for the theoretical description 
of Raman processes. 
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I. INTRODUCTION 

The spectroscopy of low-dimensional semiconductor 
structures has found in the inelastic (Raman) scattering 
of light a powerful tooli*£. By means of this technique, 
a whole sector of the electronic excitation spectrum, 
not reachable from direct light absorption because of 
Kohn theorem^, becomes the subject of research. Single- 
particle and collective states, spin- and charge excitations 
leave a trace in the Raman spectra. 

The interpretation of the results of a Raman experi- 
ment, on the other hand, is usually based on a simplified 
theoretical expression, the so called off-resonance approx- 
imation (0RA)4. The ORA states that, far away from 
resonance conditions, the Raman scattering amplitude is 
dominated by collective states, i. e., by states carrying 
on most of the strength of multipole operators. 

Experiments in quantum wells and wires under ex- 
treme resonance, however, have revealed Raman peaks 
associated with single-particle excitations (SPE) 5 . These 
peaks, absent in the ORA, are known to be related to tak- 
ing a proper account of the intermediate (virtual) states^. 
For higher excitation energies, 40-50 meV above the band 
gap, a resonant enhancement of Raman intensities for 
particular values of the incident laser energy has been 
observed^. This resonant effect is clearly not described 
by the ORA, and has been ascribed to the existence of 
incoming and outgoing resonances in the intermediate 
states. 

The relevant experimental facts concerning electronic 
Raman scattering in quantum dots can be found in Ref. 
p|, and specially in paper Q. As the incident laser en- 
ergy moves from extreme resonance to 40 meV above 
it, the observed Raman spectrum evolves from a SPE- 
dominated one to a spectrum dominated by collective 
excitations. The positions of collective excitations for 
the dots studied in Ref. have been computed in Ref. 
[i"f]| by means of ORA expressions. This is a pretty ex- 
ample in which the results are partially valid in spite of 
the fact that the physical conditions for the application 



of the ORA are not. 

In the present paper, we perform a critical examination 
of the validity of the ORA by computing Raman scatter- 
ing amplitudes for a 42-electron quantum dot. The po- 
sition of Raman peaks are then compared with the posi- 
tion of collective states, obtained from the evaluation of 
multipole operators. Let us stress that the typical dots 
used in the experiments reviewed in Q contain around 
200 electrons. We hope that our 42-electron dot captures 
the salient features of these larger systems. Notice that, 
to the best of our knowledge, Raman calculations have 
been performed only for up to 12-electron dots whithin 
a restricted Hartree-Fock schemeii. 

The plan of the paper is as follows. In the next two 
sections we include, for completeness, a summary of for- 
mulas for Raman scattering and a deduction of the ORA. 
Detailed calculations of Raman amplitudes and compar- 
ison with ORA results are presented in Sec. IIVI Con- 
cluding remarks are given in Sec. 



II. RAMAN SCATTERING BY ELECTRONIC 
EXCITATIONS 

Formulas for Raman scattering are deduced from ordi- 
nary second-order perturbation theory^. For the sake of 
completeness, we provide a summary of these formulas 
in the present section. 

The starting point is the total Hamiltonian: 



(1) 



7Yq describes the electronic subsystem in the semicon- 
ductor structure. Coulomb interactions among electrons 
are included in it, but no electron- phonon interactions 
are to be considered. This means that excitation ener- 
gies should be below the threshold for the creation of a 
longitudinal optical (Frohlich) phonon (around 30 meV 
for GaAs), which is the main electron-phonon interac- 
tion mechanism in these structures^. describes 
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the free (quantized) photon field, and Tii is the electron- 
radiation interaction Hamiltonian: 



mo is the electron mass in vacuum. Tti contains a lin- 
ear and a quadratic in A terms. In second quantization 
formalism, these terms take the form: 
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where a\ a (e\ e) are photon (electron) creation and 
annihilation operators, and e*is the light polarization vec- 
tor. 

The total wavefunction, 1^), is taken as a product of 
free photon and electronic functions. The initial state is 
made up from the ground state of the electronic subsys- 
tem and an initial photon state. On the other hand, in 



the final state the electronic subsystem is in an excited 
state, and there is a photon which energy (and wavevec- 
tor) changed from hui to hvf. The transition amplitude, 
to second order in the photon field, has contributions 

from H (A) and H^: 
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In the last equation, E refers to the energy of the elec- 
tronic system. huji„t = Ei nt — (Ei + hvi) is the difference 
between the (total) energies of the intermediate and ini- 
tial states. The sum in brackets runs over intermediate 
states giving nonzero matrix elements of 7i\ A \ In the 
experiments, the incident photon energy is close to the 
band gap energy, thus we restrict the sum to interme- 



diate states with an additional electron-hole pair, and 
will consider only the so called resonant contribution. It 
may be shown (and numerical calculations support this 
statement) that this is the main contribution to the sum. 
Taking into account Eqs. II) . and computing explicitly 
the matrix elements of photon operators, one gets: 
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In Eq. ©, Ni (Nf) represents the number of photons in the initial (scattered) beam. We will take Nf = 0. 
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Additionally, we have introduced the width, Ti„ t , of the 
level E int . P is the interband matrix element of p (con- 
sequently, in the next formulas the momentum operator 
is a dimensionless magnitude) 
by: 



and the Hf_ r are given 
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where we have introduced hole operators in the valence 
band. 

Eq. |JBJ is the main expression to be used in the fol- 
lowing sections. In Sec. IIVI we present results of Raman 
scattering in a quantum dot based on this formula. No- 
tice that the numerical evaluation of Eq. 10 requires 
obtaining reliable approximations to the electronic wave- 
functions \i), |/), \int), and the eigenvalues Ei, Ef, E int . 
In the paper, we use RPA-like wavefunctions, obtained 
with the help of methods of Nuclear Physics^, which 
are especially intended for finite systems. Details of the 
procedure can be found in Refs. 



and the spin projection is explicitly indicated. The first 
term is thus rewritten: 
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On the other hand, the sum in Eq. © has in front of 
it a big coefficient, P 2 /m . In GaAs, for example, this 
magnitude is roughly 11.35 eV—. The denominator, on 
his side, contains the difference Ei nt — Ei, which is always 
greater than the effective band gap, E gap . The sum takes 
big values when the denominator is near zero, Eint — 
Ei w hvi (resonance). To obtain the ORA, however, 
we shall assume that the photon energy is far away from 
resonance. The energy denominator is so big that we can 
neglect its variation with Ei nt . Then, we write Ei nt — 
Ei — hvi w Eg a p — hi>i, and use the completeness relation 
for the intermediate states. In this way, we get for the 
second term (the sum): 



T 2 = 



P 2 /m 
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III. THE OFF-RESONANCE APPROXIMATION 

In this section, the ORA is deduced from Eq. JBJ) in the 
limit when the incident laser energy is far from resonance. 

Let us consider the two terms in curly brackets in Eq. 
©• They are dimensionless quantities. The first term 

(A 2 ) 

comes from 7i} . It has the form of a structure function. 
Greek letters a, (3 in Eq. JfjJ denote a complete set of 
(conduction-band) one-electron states. In this section, 
we use a different notation for these states, and will write 
a f or a |, in which a labels the orbital wavefunction, 



Next, we notice that \i) and |/) do not contain hole 
states. Thus, the hole operators in Eq. (|1(J|) should con- 
tract to a Kronecker delta function. We took for the hole 
a complete basis of states labeled by the orbital index, 7, 
and its total (band) angular momentum projection, rrij. 
T2 takes the form: 



r 2 = ^{[TT] + [U] + [U] + [4T]}, (n) 
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where we have used the abbreviated notation: 
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etc. Matrix elements of factors like e ■ p depend only on 
the electron spin projection and the hole index rrij. They 
are computed according to Table [I] in which the compo- 



nents e± are defined as e± = t(^xT i£ y )/V%- Finally, 
we use completeness of the hole orbital wave functions. 
The result is as follows: 
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In the ORA, the Raman amplitude Cu is proportional 
to Ti + T 2 , where Ti and T 2 are given by Eqs. 10 and 
(|13[1 . Notice that, in this approximation, the intermedi- 
ate states play no role in the Raman amplitude. 

For simplicity, in what follows we shall consider the 
quasi-twodimensional motion of electrons. An external 
harmonic potential, with characteristic energy hu>o, con- 
fines the electron motion to a region in the plane: the 
quantum dot. 

In order to realize that only collective states contribute 
to T1+T2, let us expand the exponential e l ^~^ f ' r in Eq. 
(|13[l . The main term of this expansion, for example, leads 
to a factor: 



a, a' 
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i.e., the matrix element of the monopole operator. The 
coefficients d l ai are defined according toi£: 



d l ai = (a|plV e | 7 ); 1^0, 



(a|p 2 | 7 ); 1 = 0, 



(15) 



where p and 9 are polar coordinates in the plane. 

Final states, |/), for which D ^ is significantly different 
from zero are called collective charge monopole excita- 
tions of the initial (ground) state \i). The exact quan- 
titative meaning of this statement comes from the fact 
that D° satisfies the energy-weighted sum rulei^: 
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where [ip is the electron Fermi level. A state is (conven- 
tionally) said to be collective if its contribution to the 
sum rule is greater than 5 % of the total. Similar criteria 
can be formulated for charge multipole states (|Z| > 0)r^ 
or for spin excited states (involving or not spin reversal 
with respect to the ground state). The latter are related 
to the last three terms of Eq. 1|13[) j proportional to e*j x ef . 

We may conclude this section by saying that, in the 
ORA, only charge and spin collective states contribute 
to the Raman amplitude C/i, which is proportional to 
Ti+T 2 . 



IV. NUMERICAL RESULTS 

Below, we present results of calculations of Raman 
scattering in a quantum dot with 42 electrons. Details 
of the numerical procedure can be found elsewher o 15 ' 16 . 
We use the model parameters (effective masses, dielectric 
constant, Kohn-Luttinger parameters, etc) of Ref. |l6(. 

We shall compute the expression in curly brackets of 
Eq. ©. The first term, proportional to Si ■ e*/, shows to 
be negligible in all our computations. Thus, it will not 
be included in the following. We will also exclude the big 
factor, P 2 /mo, and will show results for the smoothed 
function 



m = E 



E 

int. 
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I is proportional to the differential cross section of the 
process. The Lorentzian in Eq. Ijl7|l is an approximation 
to the delta function coming from |C/i| 2 for large values 
of t. The phenomenological parameters, Ti nt , Tf are 
fixed to 0.5 and 0.1 meV, respectively. The sum runs 
over intermediate states with excitation energy below 30 
meV. This means states for which E int — E^ t < 30 meV, 
where E^ t is the energy of the lowest intermediate state. 
E^ t enters the definition of the qdot effective band gap: 



To approach the ORA regime, the laser excitation en- 
ergy, hvi, will range in an interval below E gap ML It will 
be shown that 30 meV below E gap is enough. 

Calculations are performed in backscattering geome- 
try, where both the incident and scattered light form an 
angle of 20° with the dot normal. 

As ORA does not provide the correct relative intensi- 
ties of the collective states, comparison between numeri- 
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cal calculations from Eq. (|17fl and the ORA expression, 
T\ + T 2l is done at a qualitative level. The position of 
collective states are obtained from the multipole opera- 
tors. When / shows to be dominated by collective states, 
we say that the ORA regime has been reached. 

We show in Fig. 2] the Raman intensity, /, in the 
charge monopole channel. This means that final elec- 
tronic states in the scattering process are monopole ex- 
citations (AL = 0) with the same total spin projection 
as the initial ground state, Computations are carried 
out in the so called polarized Raman geometry, in which 
the initial and final light polarization vectors are parallel. 
The x axis represents the Raman shift, E — Ei. The posi- 
tion of the collective state, labeled CDE (charge density 
excitation), is explicitly signalized. Peaks corresponding 
to single-particle excitations arc positioned to the left of 
the CDE. Two curves are drawn in each panel. The first 
one (shown for comparative purposes) is obtained with 
a laser excitation energy hvi = E gap + 5 meV, and the 
second one with hv\ — E gap — 30 meV. A scaling factor 
is introduced for the sake of clarity. 

We see a distinct tendency toward dominance of the 
CDE state for below band gap excitation, an indication 
that the ORA regime is being reached. This tendency 
is reinforced when the in-plane confinement, Hloq, is low- 
ered. The reason is simple. To compute /, we use in- 
termediate states in the interval \E gapi E gap + 30 meV]. 
When Hujq is decreased by a factor 77, the number of states 
in this interval increases roughly by the same factor, and 
the assumption made in the derivation of the ORA con- 
cerning completeness of the set of intermediate states is 
better fulfilled. We show in Fig. ^ curves for fkuo = 12 
and 3 meV as examples. Let us stress that the argument 
to restrict the sum over intermediate states to this en- 
ergy interval is more than technical, and is related to the 
sudden increase of Tint f° r larger excitation energies, as 
will be discussed below. 

Figs. |21 and [21 show results of calculations in other 
channels. The monopole spin-flip case, AL = 0, AS = 1, 
computed in the depolarized Raman geometry, £f£i — 0, 
is drawn in the left upper corner. The collective state, 
labeled SDE (spin density excitation), lies to the left of 
the single-particle peaks in this case. The behavior when 
hvi is varied is similar to the charge monopole situation 
discussed above. The SDE dominates the spectrum for 
below band gap excitation, and this tendency is better 
observed for hwo = 3 meV. 

In the charge quadrupole channel, AL = 2, AS = 0, 
and polarized Raman geometry, right lower panel, the 
dominance of the CDE is observed only when fkvo is low- 
ered to 3 meV, accompanied by a strong decrease of the 
intensity. 

Dipole final states, AL = ±1, however, show a dif- 
ferent behavior in which collective states are never the 
dominant peaks in Raman scattering. The charge dipole 
channel, AL = 1, AS = 0, computed in the polarized 
geometry, is shown in the right upper panel. Notice that 
the CDE is not a significant peak even with above band 



gap excitation. In the spin-flip case, left lower corner, 
the SDE state is always a low peak. Notice, in addition, 
that in the latter case the intensity experiences a strong 
decrease for below band gap excitation. 



V. DISCUSSION 

We have computed Raman intensities from a 42- 
electron qdot in different angular momentum and spin 
channels for below band gap excitation to show how and 
to what extent the ORA regime is reached. 

It was demonstrated that for (charge and spin) 
monopole (AL = 0) final states, collective excitations 
rapidly become the dominant excited modes and, already 
for laser energy 30 meV below the effective band gap, are 
the main peaks in the Raman spectra. 

In any other channel, however, the ORA may or may 
not play a relevant role. We observed, for example, the 
dominance of the CDE state in the charge quadrupole 
channel when Hloq is lowered, and thus the density of 
states in the dot is increased. But we observed no range of 
hvi in which collective modes become the leading dipole 
peaks. 

The conclusion is that the ORA works very well for 
monopole states, but only conditionally for quadrupole 
ones, and very badly for dipole states. The number of 
electrons in the dot or the strength of the external con- 
fining potential may dictate the way in which the ORA 
regime is reached by varying the density of energy levels 
in the dot. Experiments in quantum dots for below band 
gap excitation could confirm this conclusion. 

More generally speaking, we suggest caution when an 
experimental Raman peak is to be interpreted as a col- 
lective excitation, on the basis of ORA arguments, in 
other systems like quantum wires or wells. Theoretical 
calculations should support this kind of assignements. 

Collective states could dominate a Raman spectrum 
in a regime very different to that one leading to the 
ORA. Let us recall again the results presented in pa- 
per |9j for laser excitation energies 40 - 50 meV above 
the effective band gap. In this energy interval there are 
plenty of intermediate states, thus the experiment deals 
with a resonant regime. A correct theoretical description 
should, probably, take account of the dependence of Ti nt 
on Ei nt i^ On qualitative grounds, one can expect a sud- 
den increase of Ti nt when Ei nt is over 40 meV above the 
effective band gap. There could be, however, special in- 
termediate states (which can be thought of as "plasmon 
plus exciton" states) for which Ti nt may preserve rela- 
tively small values. The physics of Raman scattering in 
this regime is based, in our opinion, on the existence of 
these intermediate state resonances, and has no relation 
with the ORA. This interesting problem deserves more 
theoretical and experimental work on it. 

A second example of collective states dominating Ra- 
man scattering, and in which the ORA has no relevance, 
is given in Ref. 19], where the peak related to the elec- 
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tronic cyclotron resonance in a quantum well undergoes 
an abrupt increase when it moves through the position 
of a luminescence peak (this is achieved by varying hvi). 
The theoretical description, we guess, should include the 
enhancement of Raman intensity due to the background 
luminescence (The Nf variable, which was set to zero 
in Eq. ©) and a consistent calculation of luminescence 
(including master equations for the population of the in- 
termediate states). This situation is, again, far from the 



simple assumptions leading to the ORA. 

Our final conclusion is, thus, that there are still many 
interesting experimental and theoretical aspects of Ra- 
man scattering in semiconductor quantum structures 
waiting for elucidation. Simple recipes such as the ORA 
give a first insight into the physics, but more elaborated 
theoretical constructions are needed for a complete un- 
derstanding of the process. 
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TABLE I: The matrix elements {a\e •p\mj)/i 
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FIG. 1: Raman intensities for charge monopole final states and polarized geometry: (a) huo — 12 meV, (b) hwo = 3 meV. Two 
curves are shown corresponding to hfi = E gap + 5 meV (solid line), and hi>i = E gap — 30 meV (dashed line). This convention 
is used in the next two figures. 
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FIG. 2: Raman intensities in different spin and angular momentum channels, and for hujo —12 meV. 
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